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Abstract. 

The mass distribution of dark matter halos is a sensitive probe of primordial 
non-Gaussianity (NG). We derive an analytical formula of the halo mass function by 
perturbatively computing the excursion set path-integrals for a non-Gaussian density 
field with no n- vanishing skewness, /jyx- We assume a stochastic barrier model which 
captures the main features of the ellipsoidal collapse of halos. Contrary to previous 
results based on extensions of the Press-Schechter formalism to NG initial conditions, 
we find that the non-spherical collapse of halos directly alter the signature of primordial 
NG. This points toward a potential degeneracy between the effect of primordial non- 
Gaussianity and that of non-linear halo collapse. The inferred mass function is found 
to be in remarkable agreement with N-body simulations of NG local type. Deviations 
are well within numerical uncertainties for all values of —80 < f l § c L < 300 in the range 
of validity of the perturbative calculation. Moreover, the comparison with simulation 
results suggests that for f NL > 150 or f NL < —50 the non-linear collapse of halos, 
as described by our barrier model, strongly deviates from that of Gaussian initial 
conditions. This is not surprising since the effect of non-linear gravitational processes 
may be altered by initially large NG. Hence, in the lack of prior theoretical knowledge, 
halo collapse model parameters should be included in statistical halo mass function 
data analysis which aim to constrain the signature of primordial NG. 
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1 Introduction 

The statistics of the primordial density fluctuation field carries unique information on 
the physics of the very early Universe. Standard slow-roll inflation predicts a nearly 
scale invariant spectrum of Gaussian, adiabatic density perturbations [1-7]. Alterna- 
tively, high-energy physics inspired scenarios may leave large non-Gaussian signatures 
on the initial matter density distribution [8-19]. Hence, the detection of primordial 
non-Gaussianity (NG) has the potential to disclose quantum physical processes at 
energies which are out of reach of laboratory experiments. In recent years this has mo- 
tivated a revived interest on testing the Gaussian hypothesis through measurements 
of the high-order angular correlation functions of the Cosmic Microwave Background 
(CMB) anisotropies [20-23]. 

The quest for primordial NG has primarily focused on the amplitude and shape 
of the CMB bispectrum statistics [24-27]. The next generation of cosmic structure 
surveys may also provide complementary constraints, since the spatial distribution of 
galaxies in the universe carries an imprint of the statistics of the initial matter density 
fluctuations. As shown in several studies, measurements of the galaxy bispetrum [28- 
32] and the abundance of galaxy clusters [33-35] are sensitive NG probes. However, 



- 1 - 



an accurate modeling of the non-linear regime of gravitational collapse is necessary for 
such tests to be successful. This is because the late time non-linear evolution of density 
perturbations induces strong non-Gaussian features in the spatial distribution of cos- 
mic structures, thus potentially obscuring the imprint of primordial non-Gaussianity. 
Accounting for these effects is the premise to testing the statistics of the initial density 
field. 

In this context the determination of the mass distribution of dark matter halos 
is of primary importance. Halos are gravitationally bound objects resulting from the 
non-linear collapse of dark matter perturbations. It is inside these virialized struc- 
tures that cooling baryonic gas falls in to form the stars and galaxies that we observe 
today. Numerical N-body simulations have been extensively used to study the non- 
linear clustering of matter in the standard cosmological scenario with Gaussian initial 
conditions. On the other hand, simulations with non-Gaussian initial conditions have 
become available only in the past few years (e.g., see refs. [36-39]). This has driven 
a major effort toward formulating a mathematical model description of the NG halo 
mass distribution. A number of studies have approached this problem by extending 
the original Press-Schechther (PS) derivation to NG initial conditions. The seminal 
work by Press and Schechter [40] has provided the first formal derivation of the halo 
abundance in the case of an initial Gaussian density field. The basic idea is that halos 
form in regions of the linear density field smoothed on a given scale R (associated to a 
mass M), whose density lies above a critical threshold of collapse (e.g. 5 C the linearly 
extrapolated spherical top-hat density perturbation at the time of collapse [41]). Using 
such a prescription, the fraction of mass in halos with mass > M is given by 

POO 

F(> M) — I U(5,M)d5, (1.1) 
Js c 

where II (5, M) is the probability distribution of the smoothed linear density field. 
Then, the number density of halos with mass in the range [M, M + dM] is obtained 
by equating 

dM MdM' 1 ' ' 

with p being the mean background matter density. 

The computation can be extended to non-Gaussian initial conditions by specifying 
the form of 11(5, M). For example, Matarrese, Verde and Jimenez [42] have written 
11(5, M) as a path-integral of an infinite series of cumulants of the density field 5. This 
allowed them to derive an approximate formula for dn/dM in the limit of small NG 
up to leading order in the skewness parameter, /jvl- Similarly, in [43] the authors have 
approximated II(<5, M) as an Edgeworth expansion in the cumulants and truncated the 
series to infer an expression of the mass function which explicitely depends on f^L (for 
an extension to the ellipsoidal collapse threshold see ref. [44]). Overall these studies 
have provided indications of how primordial NG affects the halo mass distribution. 
However, by relying on the PS approach they have failed to address a number of 
important issues that are the main limitation of the PS formalism. As an example in 
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Press-Schechter the normalization factor of the mass function is related to the "cloud- 
in-cloud" problem. Namely, the fact that a collapsed region of mass M contributes 
to the halo mass counting only if it is not embedded in a collapsed region of greater 
mass M' > M. The PS formalism does not address this issue, thus miscounting the 
fraction of mass in halos. Moreover, even including an ad-hoc correction to recover 
the correct normalization, the PS formula still fails to reproduce results from N-body 
simulations. Because of this, analyses based on the PS approach have mainly focused 
on predicting the ratio of the non-Gaussian to Gaussian Press-Schechter mass function 
and normalized the result with a fitting function to Gaussian N-body simulation data. 
A recent comparison with non-Gaussian simulations indicates that these analytical 
formulae are in relatively good agreement with numerical results provided an ad-hoc 
rescaling of the spherical collapse threshold is assumed [37]. 

Ideally, one would like to perform an ab initio calculation to consistently account 
for the cloud-in-cloud problem as well as the non-spherical collapse of halos for any 
type of initial conditions. The excursion set theory introduced by Bond et al. [45] 
provides a powerful mathematical framework capable of addressing these issues and 
allowing for a self-consistent derivation of the mass function from a limited set of 
initial assumptions. The formalism frames the original Press-Schechter idea in the 
context of stochastic calculus. However, the difficulty in obtaining analytical solutions 
for realistic description of halos has been the main limitation of the theory (for a 
review see ref. [46]). In fact, in the case of a Gaussian field an analytic treatment is 
possible only if one assume an unphysical halo mass definition. The computation for 
a realistic case requires to solve numerically the stochastic model equations through 
time consuming Monte Carlo simulations (see e.g. refs. [45, 47]). 

In a series of papers Maggiore & Riotto [48-50] have introduced a path-integral 
description of the excursion set which allows for a perturbative calculation of the mass 
function in the case of physical mass definitions. Non-Gaussian initial conditions can 
be easily implemented in this novel formulation as well as more realistic conditions for 
the gravitational collapse of halos. This has provided a major breakthrough towards an 
accurate physical model description of the mass function. For instance, in [51, 52] we 
have implemented such an approach with an effective stochastic barrier model which 
accounts for the main characteristics of the ellipsoidal collapse of halos. This has 
allowed us to derived an analytic formula which shows an unprecedented agreement 
with the mass function from Gaussian N-body simulations. 

In this paper we extend the study of the non-spherical collapse mass function to 
the case of primordial non-Gaussianity. Our goal is to derive a better understanding of 
how the halo collapse alters the signature of primordial NG and confront the theoretical 
model prediction against non-Gaussian N-body simulation results. The paper is orga- 
nized as follows: in section 2 we introduce the excursion set theory, the path-integral 
formalism and discuss diffusive drifting barrier model of halo collapse. In section 3 we 
derive the non-Gaussian mass function and shown its application to local and equilat- 
eral non-Gaussianity in section 4. We discuss the imprints on the halo mass function 
due to primordial NG and the halo collapse model in section 5, while in section 6 we 
present the results of the comparison to non-Gaussian N-body simulations. Finally we 
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discuss the conclusions in section 7. 



2 Excursion Set Theory 

2.1 Correlated Random Walks and Halo Mass 

Let us consider a non-Gaussian density field 5(x) smoothed on a scale R: 



<J(x, R) = y d 3 k W(k, R)~5(k)e-* k x (2.i; 



where W(k,R) is the Fourier transform of the filter function in real space, W(x,R). 
The latter defines the mass inside a region of radius R as M = pV(R), with p being 
the mean background density and V(R) = J d 3 x W(x, R) is the filtered volume. Then, 
the variance of the smoothed density field reads as 

S = (5 2 (R)} = a 2 (R[M\) = ^ J dk k 2 P(k)W 2 [k, R(M)\, (2.2) 

where P(k) is the linear matter power spectrum at redshift z = 0. Given the one-to- 
one relation between smoothing radius, mass and variance, we can rewrite the halo 
mass function eq. (1.2) as 

dn _ t( . p rflogq- 1 

dM~ na) M 2 rflogM' {Z - 6) 

where f(a) = 2a 2 J r (a 2 ) is the so called "multiplicity function" and J-{S) = dF/dS 
is the derivative of the fraction of mass in halos with mass > M[5]. The goal of the 
excursion set theory is to calculate f(a) from F(S) obtained as the first up-crossing 
distribution of the random walks 5(x, R) as function of the smoothing radius R[S] (or 
equivalently S, which plays the role of pseudo-time). Trajectories starts at S = (i.e. 
large radii limit R — > oo) with 5q = 0. The fraction of mass in halos is then inferred from 
counting the rate at which trajectories hit for the first time an absorbing barrier. The 
mass function at high redshifts can be inferred by simply rescaling all linear quantities 
by the growth factor predicted by the cosmological model under consideration. 

The nature of the random walks depends on the statistics of the initial density 
field as well as the form of the filter function. In the case of a sharp-k filter, W(k, R) = 
0(1 — kR), the random walks are Markovian with (S(S)S(S')) C = min(5', S'). 1 In such 
a case the first-crossing distribution is obtained by solving a simple Fokker-Planck 
diffusion problem. However, the sharp-k filtering does not correspond to a physical 
halo mass, since the filtered volume V(R) is undefined (see e.g. discussion in [48]). 
The only unambigously defined mass is associated with a top-hat filter in real space, 
W(x, R) = 6(R — x) (sharp-x filter) whose Fourier transform reads as 

W(k, R) = -^—\sm(kR) - (kR) cos(kR)]. (2.4) 
(kR) 6 



1 We use the underscore c to indicated connected correlators. 
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In such a case the evolution of the random walks acquires pseudo-time correlations for 
which to our knowledge an analytical solution of the first-crossing distribution does 
not exist. Maggiore & Riotto [48] have shown that such correlations can be treated 
perturbatively around the Markovian solution. In particular, for a sharp-x filter the 
two point connected correlator can be written as (6(S)5(S')) C = mm(S, S') + A(S, S'), 
where A(S,S') ~ kS/S'(S' — S) < 1 (with S < S') and k is a small parameter of 
nearly constant amplitude, whose value depends on the underlying cosmology. For a 
vanilla ACDM model k ~ 0.475. Hence, the non-Markovian first-crossing distribution 
can be inferred perturbative expansion in k. 

2.2 Stochastic Barrier and Ellipsoidal Collapse 
2.2.1 State of art of the stochastic barrier model 

The absorbing barrier is the other crucial ingredient of the excursion set formalism 
which encodes information on the non-linear collapse of halos. In the spherical col- 
lapse model, the barrier is a constant threshold, B = 5 C . However, this model is a 
too simple description of the collapse of dark matter halos in the universe. The grav- 
itational collapse can be highly non-spherical. As already shown in the seminal work 
by Doroshkevich [60] perturbations in a Gaussian density field are triaxial and subject 
to external tidal shear. The implications of this are twofold. First, the collapse of a 
homogeneous ellipsoid [61] provides a better description of the system. Secondly, the 
stochastic nature of the underlying density field implies that the parameters character- 
izing the ellipsoid (ellipticity and prolatness) are random variables with a well defined 
probability distribution that can be derived from that of the density field. Conse- 
quently, the condition of collapse is no longer specified by a unique constant value, 
but it is a random variable itself [53, 54]. In the excursion set framework this implies 
treating the barrier as a stochastic variable which performs a random walk whose char- 
acteristics are given by the moments of the probability distribution resulting from the 
ellipsoidal collapse model considered. Unfortunatly the full distribution of the barrier 
remains unset. Most of the work dedicated to the ellipsoidal collapse threshold have 
limited their analysis to the average behavior of the collapse threshold as function of 
the halo mass. As an example Sheth, Mo & Tormen found that the average ellipsoidal 
collapse barrier evolves as (B(S)) = 8 C (1 + /3S" 7 ). This behavior is consistent with the 
physical expectation that collpase of small mass halos differ from spherical symmetry. 
In fact, at large S (small mass) the shear field is higher than at small S (large mass) 
thus implying that small mass halos need to reach a higher density of collapse than 
that predicted by the spherical model. However, knowledge of the scatter around this 
average encodes additional information on the non-linear collapse process. To date the 
variance of the stochastic barrier has been studied only qualitatively [54]. In [62] it has 
been shown that the dynamics of the collapse depends not only on the local properties 
of the shear field but also on interaction between ellipsoidal perturbations and its large 
scales environment. One of the key parameter to quantify this effect is the large scale 
environment density. Low mass haloes are more affected by this large scale environ- 
ment and one of the consequence is to increase the scatter. This study was done for a 
Gaussian field but it already suggests that the imprint of the primordial density field 
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is underline at hight mass halo while the dynamics and the cosmology dependences 
induce by the environment dynamics is underline in low mass haloes. 

However we must not draw any conclusion on the value of these parameters with- 
out a proper account of the second order cloud-in-cloud problem. In the excursion set, 
we assume that we can consider random walks around the same coordinate holding on 
the ergodic principle. This hypothesis is approximative and it is possible to measure 
the induced discrepancy on the scatter by doing the analysis of [54]: comparing the 
scatter for center of mass walk and random particle of the simulation. Nevertheless 
this will require a dedicated study that we would like to perform in the future including 
non-Gaussian statistic of the primordial field. In the absence of prescription we focus 
a simpler model with captures the main feature of the fuzzy barrier. 

2.2.2 Modelisation of the Diffusing Drifting Barrier 

As we already introduce in [51, 52] we use a stochastic barrier with Gaussian diffusion 
and linear drifting average, which we refer as the "Diffusive Drifting Barrier" model 
(DDB). In this model the barrier performs a Markovian random walk starting at B = 
5 C with (B(S)) C = 8 C + (3 S and (B(S)B(S')) C = Db min( f S', S'). Here (3 parametrizes 
the rate at which the collapse threshold on average deviates from the spherical collapse 
prediction and Db measures the amplitude of a constant scatter around the average 
collapse threshold at a given mass scale. The linear behaviour of the average of the 
threshold approximates the mean value of the Sheth Mo and Tormen ellispoidal barrier 
in the mass of interest. 

Therefore, differently from the linear density field 5, the non- Markovian part of 
the barrier 2-point connected correlator identically vanishes, Ab(S, S') = 0. One may 
wonder if a better modeling of the stochastic barrier should account for the correlations 
induced by the filtering process. However, as already pointed out in [52], there is no 
reason as to why the barrier should have the same filtering of the linear density field, 
since the two procedures have very different physical meanings. The latter is associated 
to the halo mass definition, while the former specifies the correlation between the 
condition of collapse at different scales. To date there is no study that has looked 
at such correlation in the context of the ellipsoidal collapse model. In the lack of 
information, the simplest approximation is to assume that the collapse condition at 
scale S is independent of that at S', i.e. Ab(S,S') = 0. This suggests that in order 
to improve the modeling of the stochastic barrier one should completely solve the 
ellipsoidal collapse model through a detailed numerical analysis from which to infer 
moments and correlators of the fuzzy barrier. This will also allow to predict the 
redshift and cosmology dependence of (3 and D B - 

The stochastic system can be simplified into one-dimensional random walks by 
introducing the variable Y = B — 5, with initial condition Yq = 5 C and absorbing 
boundary at Y c = 0. Non-Gaussian initial conditions can be easily included through 
the higher-order connected correlation functions of the smoothed linear density field. 
For instance, let us consider non-Gaussianity generated by the three-point correlation 
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function (skewness). In such a case the random walks are completely specified by the 
non- vanishing connected correlators 

<y{S)) c = s e + ps, (2.5) 

{Y(S)Y(S')) C = (1 + D B )mm{S, S') + A(S, S'), (2.6) 
(Y(S)Y(S')Y(S")) C = -{5(S)5(S')5(S")) C} (2.7) 

where (5(S)5(S')5(S")) C is the three point correlation function of the smoothed linear 
density field. The goal of the excursion set is to compute H(Y , Y, S), the probability 
distribution of the random walks which start at Y and reach Y at time S obeying 
eq. (2.5), (2.6) and (2.7), with U(Y ,Y C , S) = 0. Then, the first-crossing distribution 
can be derived as 

= -lL [°° dY n ( y o, Y, S). (2.8) 



dS 



y, 



2.3 Path-Integral Formulation 

The computation of n(y ; Y, S) can be performed using the path-integral method in- 
troduce in [48]. Let us consider a discretized interval [0,S] in steps AS" = e, such that 
Sk = ke with k = 1, ...,n. The probability distribution of the discrete random walks 
that never cross the barrier reads as 

n e (F , Y n , S n ) = dY!... dY n ^W(Y , .., Y n , S n ), (2.9) 
Jy c Jy c 

where 

W(Y , .., Y n , S n ) = (5d(^(S'i) - Y 1 )...5 D {Y{S n ) - Y n )), (2.10) 

is the probability density distribution. Using the Fourier transform of the Dirac- 
function we can rewrite eq. (2.10) as 

W(Y ,..,Y n ,S n ) = J DAe'^ A ' y '(e- ! ^ A ' Y(Sl) ), (2.11) 

with J T>\ = ^-•••^f • As clearly pointed out in [50], the exponential term 
inside the average is nothing else than the exponential of the partition function, 
e z = (e -4 ^^ 1 A * y ( 5 *)). For generic non-Gaussian random walks this reads as 

oo / ,\„ n n 

z = E E - E ^-KiY^.YiX (2.12) 

p=l ii=l i p =l 

where (Y il ...Y i ) c is the p-point connected correlator. Using eq. (2.5), (2.6) and (2.7) 
and substituting in eq. (2.10) we have 

/•oo /*oo /' 

U e {Y ,Y n ,S n ) = / dY x ... / dY n ^ / V\ e^M^-nH^A^. x 
Jy c Jy c J 

x e~l Ei i AiA J A «i e ^-r- s ijfc A i A J A fc< y i Y j y fc)c ; (2.13) 
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with Aij = e(l+Z?B)min(z, j). Since Ay = A(jSj, Sj) < 1 and assuming small departure 
from Gaussianity, we can expand the exponentials in eq. (2.13) to first order in A^ and 
{YjYjY k ) c respectively. Then using the relation \ie~ tX ^ Bi ~ Y ^ = —id / dYie~ tT,kXk ^ Bk ~ Yk ^ 
we can write 

n £ (F , Y n , S n ) « W (Y , Y n , S n ) + Ul(Y , Y n , S n ) + U% G (Y 0l Y n , S n ), (2.14) 

where 

poo poo 

W (Y , Y n , S n ) = dY 1 ... / dY n -iW (Y , Y n , S n ), (2.15) 
Jy c Jy c 

is the Markovian probability distribution with density 

W (Y ,...,Y n ,S n ) = J v\e~^ kXk[Bk - Yk] -^ x ^ A >i, (2.16) 
which obeys the Chapman-Kolmogorov equation (see appendix A in ref. [52]), 

1 poo poo 

ILXXo, Y n , S n ) = -J2 dY ±- / dY -^ AijdtdjWoiYo, .., Y n , S n ), (2.17) 
is the first order correction in k and 

1 POO POO 

IL% G (Y , Y n , S n ) = -- J2(^ 3 5 k ) c / dY 1 ... / dY n ^ d&dkWoiYo, .., Y n , S n ) 

ijk ^ Yc 

(2.18) 

is the non-Gaussian term [50]. Notice that due to the separate expansions in A^ and 
(YjYjYk) c we have mixing terms. However, as we will discuss later, these are negligible 
to first order in k and /jvl- By evaluating these integrals in the continuous limit 
(e — > 0) and computing the corresponding first-crossing distribution we obtain a mass 
function that consists of three independent terms 

f^) = f {a) + f 1 (a) + f NG {a), (2.19) 

where /o(c) is the Markovian part of the Gaussian mass function, fi(a) is the non- 
Markovian correction to the Markovian solution to first order in k due to the filtering 
process of the linear density field and /jvg(c) is the non-Gaussian contribution. Before 
presenting the computation of /atg(c)) it is worth reviewing the main features of the 
Gaussian mass function derived in [51, 52]. 

2.4 Gaussian Halo Mass Function 

The Markovian mass function for a diffusive drifting barrier reads as 

fo(cr) = 6 -lM e -^ 5 ^ 2 . (2.20) 
a V 7r 

where a = 1/(1 + Db)- The non-Markovian correction due to the filter function is 
given as an expasion to second order in f3 

hi?) = /!7=W + + (2-2i) 
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with 





fm—m/\ 
71,0(2) 




(2.22) 



(2.24) 



(2.23) 



where Yq = 5 C and T(0, z) is the incomplete Gamma function. This expansion is 
justified by the fact that from the ellipsoidal collapse model we expect < 1 and in 
[51, 52] we have shown that terms 0(> /3 3 ) are negligible 2 , 

The collapse model affects the Gaussian mass function, /g(ct) = /o(c) + f\{o), 
in a non-trivial manner (for an exhaustive discussion we refer the reader to section 
III in ref. [52]). In figure 1 we plot /g(c) f° r different values of /3 and Db over the 
interval corresponding to the mass range probed by current NG N-body simulations 
(10 13 < M[h _1 M ] < 10 15 ). The spherical collapse scenario studied in [48] corresponds 
to (3 = and Db = 0. We can see that the effect of the average drift of the barrier 
(see curve with (3 = 0.3 and D B = 0) is to suppress the mass function at all masses 
with respect to the spherical collapse case. Since the average threshold of collapse 
drifts toward larger values as function of a, it becomes less likely for random walks to 
first-up crossing the barrier, thus decreasing the amplitude of /g(°") more effectively in 
the low-mass end than at larger masses. On the other hand the effect of the stochastic 
diffusion (see curve with (3 = and D B = 0.3) is to facilitate the barrier crossing, thus 
effectively lowering the threshold of collapse. This is consistent with the presence of 
the factor a — 1/(1 + Db) in the exponential cut-off of the mass function, which leads 
to a larger amplitude of the mass function especially in the high-mass end. In the 
case of a diffusing drifting barrier (see curve with j3 = 0.3 and D B = 0.3) these two 
effects compete to produce a mass function which is tilted with respect to the spherical 
collapse prediction. 

The inclusion of the non-Markovian correction f\ (a) due to the filtering procedure 
has the effect of decreasing the amplitude of the mass function. This is because in the 
presence of pseudo-time correlations the random walks are characterized by smoother 
trajectories with less frequent jumps. This has the effect of disfavoring the first-up 
crossing of the barrier, thus lowering the mass function. This is consistent with the 
overall minus sign in eq. (2.22) and eq. (2.23) which are the leading and next-to- 
leading non-Markovian corrections in the mass range of interest. Here, it is also worth 

2 In f^g^icr) we had previously miss a contribution. However this term is negligeable for mass 



less than 10 16 M SOTl in [51, 52] 
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0.2 0.4 

log 10 (CT" 1 ) 



Figure 1. Gaussian excursion set halo mass function with sharp-x filter for the spherical 
collapse barrier (black solid line: /3 = 0, Db = 0), a spherical diffusive barrier (green dot 
line: (3 = 0, Db = 0.3), rigid moving barrier (blue short dash line: /3 = 0.3, Db = 0) and 
diffusive drifting barrier (red long dash line: /3 = 0.3, Db = 0.3). 

noticing that these expressions carry higher powers of a compared to the Markovian 
term fo{&)- This is because the effect of pseudo-time correlations on the mass function 
is diluted in the presence of diffusion (since diffusion always favourates barrier crossing), 
thus leading to a faster convergence of the expansion in k. As shown in [52], the 
mass function with non-Markovian corrections to first order in k is consistent within 
numerical errors with the exact numerical solution inferred from Monte Carlo generated 
random walks. 

The study presented in [51, 52] has shown that implementing a realistic modeling 
of the non-spherical collapse of halos in the excursion set formalism is key to inferring a 
physical description of the halo mass function which is capable of reproducing results 
from N-body to numerical accuracy. Here, we will show that this is the case also 
for non-Gaussian initial conditions, since the halo collapse model directly alters the 
non-Gaussian signal. Previous works on the non-Gaussian halo mass function have 
simply reabsorbed the effects of the ellipsoidal collapse of halos in the overall Gaussian 
normalization function, for instance using the Sheth-Tormen formula [55]. But in doing 
so they have implicitly assumed that the collapse is independent of the statistics of 
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the primordial density field and thus does not affect the non-Gaussian contribution to 
the halo mass function. However, we have seen that the barrier model affects both the 
Markovian and non-Markovian part of the halo mass function. Since non-Gaussianity 
appears as a non-Markovian correction to the Gaussian term, it is reasonable to expect 
that the stochastic barrier directly affects the dependence of the halo mass function 
on primordial non-Gaussianity. We will show this in detail in the next sections. 



3 Non-Gaussian Halo Mass Function 

In order to compute eq. (2.18) we need to first address the scale dependence of the 
three-point correlator of the smoothed linear non-Gaussian density field. Following 
[50] we expand (5i5j5k) c in a triple Taylor series around Si = Sj = Sk = S n 3 : 

(5(s t )5(s 3 )5(s k )) c = J2 , , r , (s n -s i Y(s n -s j y(s n -s k yc$«' r \s n ), (3.1) 

p,q,r=0 P-1- r - 



where 

, s rlP dl ,V 

(3.2) 



G { r r \s n ) = ^^^(SMSMSk)), 



dSi dSj dS r k s . =Sj=Sk=Sn 

In the large mass limit (S n —> 0), where we expect that the primordial NG signature 
to be less altered by non-linear effects, the leading contribution is given by the term 
(S 3 (S n )), corresponding to p + q + r = 0. The next-to- leading order is given by 
p + q + r = 1 corresponding to three terms proportional to derivative of (5(Si)5 2 (S n )) 
and so on, thus 

^ng(Yo, Y n , Sn) — n% G (Yo, Y n , S n ) + (Yq, Y n , S n ) + ... (3-3) 

In [57], the authors have computed the excursion set non-Gaussian mass function for 
a generic rigidly moving barrier modeled as a Taylor expansion in S n up to leading 
order. Here, we provide an exact computation for the stochastic barrier with linearly 
drifting average up to next-to-leading order. 

3.1 Leading Term 

The zero order term p = q = r = reads as 

n /"CO /*oo 

U% L G (Y ,Y n ,S n ) = --(5 3 (S n )) / / dY^didjdkWoiYo^.^Sn), 

(3.4) 

We can decompose the sum as 

n n—1 n—1 n—1 

didA = a> + 3£a£a + 3X; + £ a^a*, (3.5) 

i,j,k=l i=l i)i=l i,j>k=l 



3 In [56] the authors compute the path-integral through a saddle-point technique and identify scale 
dependent small parameters which allow them to infer a mass function which is valid over a larger 
range of masses. 
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the integrals associated to these terms can be computed using the fact that 



_d_ 



™ 1 poo poo 

I$(Y , Y n , S n ) = -J2 dY x ... / dY^diWoiYo, .., Y n , S n ) 



& 1 



n-l 



dY? 



I^(Y ,Y n ,S n )=J2 I dY,... I dY r ^ 1 d i d j W (Y (h ..,Y n ,S n ), 



y c 



(3.6) 
(3.7) 



and 



d 3 
my 



,L poo poo 

n^(F , Y n , S n ) = -J2 / dY i- / dY n-i d t d 3 d k W {Y , .., Y n , S n ). (3.8) 

i,j,k=i ^ Yc 

The first-crossing distribution is given by eq. (2.8) and one can easily show that the 
integrals in dY n (in the continous limit) of the first three terms in eq. (3.5) vanish. 
Thus, the only non-vanishing contribution to the first- crossing distribution is given by 
the term in eq. (3.8) such that: 4 



(S 3 (S)) 



dY? 



where 



poo 

I L = / IL {Y ,Y c ,Y,S)dY. 
Jy c 

This integral can be computed exactly using the Markovian solution [52] 



U (Y ,Y C ,Y,S) 



2-kS 



o aP(Y-Y -PS/2) 



e -^siY-Y f _ e -^(2Y c -Y-Y f 



obtaining 



h = -{ Erfc 



( Y c - Y - PS 

\y/2S(l + D B ) 



- e p 1+d b Erfc 



/ Y -Y C -PS 

\y/2S(l+D B ) 



(3.9) 
(3.10) 

(3.11) 
(3.12) 



Substituing this expression in eq. (3.9) and using the fact that f(a) = 2<j 2 J r (cr) we 
finally obtain 

R a(Y + (3a 2 ^ 2 

fNaW 



2 - 

7T 



2a 2 



a 3 / 2 a 



6 



-F 4 - 2-Y* - 1 + 3-/3F 



a' 



a- 1 



a 



+3a(3Y + a 2 (3 3 a 2 Y + 3a 2 ((3Y ) 2 + 13a/?V] + 
alner Vo l 



(3.13) 



+ e^ 2a/3y °Erfc 



2a 2 



(Y - (3a 2 ) 



2a s a 



(3 3 a 3 45 3 (a) 



dS 3 (a) 
dlna 



4 Following the same derivation it is easy to obtain the non-Gaussian contribution to the halo mass 
function from a non- vanishing 4-point correlation function by expanding the path- integral eq. (2.11) 



to first order in (SiSjSkSk)- Then, to leading order in Si = Sj = Sk = Si = S n the first crossing 

>(s))^r/~dyn (y ,y,5)' . 



distribution is — jf^jj 
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where we have introduced the reduced cumulant 63(17) = (5 3 (S)) /a 4 . For the spherical 
collapse barrier, i.e. (3 = and a = 1 (Db — 0), the above expression reduces to 
Loverde et al. [43] formula as already noticed in [50]. On the other hand, it is worth 
noticing the presence of additional mixing terms proportional to powers of (3 which 
couple the effects of the non-spherical halo collapse to primordial non-Gaussianity. 

3.2 Next-to-Leading Term 

The next term in the expansion eq. (3.1) correspond to p + q + r = 1 (i.e. p = 1 and 
q = r = plus permutations) 



1 n " POO POO 

Tl% N G L (Y , Y n , S n ) = 3x -G$ m (S n ) J2( S n " Z / dYl - / W 

i=l j,k=l^ Yc ^ Yc 



%W , 
(3.14) 



where 5 



G 



(1,0,0) 



{S n ) = 



(3.15) 



Si — Sn 



Again, decomposing the sum and using £•;=! -> d 2 /dY c 2 , we obtain 



poo 1 _ _ ffi r /"°° z" 00 

/ dY n n% N J J (Y ,Y n ,S n ) = -G^ 0fi \S n )Y / (S n -S i )— / d^... / dY n cW, 

(3.16) 

The multiple integral in the above expression can be computed by part and using the 
fact that Wo obeys the Chapman-Kolmogorv equation we can rewrite the integrand as 

Wo(^0) ••) Yi-i, Yi — Y c , Yi + i, Y n , S n ) = Wo(Yo, .., Y c , Si)Wa(Y c , Y n , S n — Si). 

(3.17) 

After taking the continous limit, the first-crossing distribution at next-to-leading order 
reads as 



3~ ng(S) 



ld_ 

2dS 



G 



(1,0,0) 



(S) 



d 2 i 



NL 



dY} 



where 



' NL 



pS poo 

/ dS'(S-S')U f (Y ,Y c ,S') dYU f (Y c ,Y,S-S'), 
Jo Jy c 



(3.18) 



(3.19) 



with IIq(Yo) Y c , S) and Hq(Y c , Y,S — S') are the finite pseudo-time corrections at the 
barrier location (see [52]): 

U f (Y , Y c , S) = ~^(Y C - Y o )e-*0ro-y^P3f (3-20) 



n f o(Y c ,Y,S) 



^s 3 / 2 

a 

^S 3 ' 2 



(Y c - Yo )e-^ Y °- Y ^ 
(Y-Y c )e-^ Y ~ Y ^\ 



(3.21) 



5 We have G^"°'°\S n ) = G^' 1 '"' (S n ) = G^'' 1 '' 1 ' (S n ) since the 3-point correlation function (Si5j5k) 
is symmetric under permutation of the indeces, a property that results of the fact that its dependence 
on Si, Sj and Sk is given by the product of three identical filter functions computed at Rj and 
Rk respectively. 



-» (0,1,0), 



■>(o,o,i). 
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eq. (3.19) can be computed analytically using the fact that 



e 2s ie 2 (S - Sl) i 2n _± (a+b) 2 
o sl^S-Si)*/* 6 V S 



we obtain 



I NL = a((3a 2 +Y c -Y )e 2a ^ Yc - Yo) Edc 



2a 2 



(Y -Y c - (3a' 



l 2 r~ 
+\j -y/aae 



(3.22) 



(Y -Y c + ^a 2 ) 2 
2a 2 

(3.23) 



Introducing Us(S) = 4 Gg (5 1 ) we finally obtain the multiplicity function to 



next-leading order 



J NG ( ) 



1 2a 2 



c?l 2 a 



1 + — Y 2 + 15af3 2 a 2 + Aa(3Y 



a 



+ ^M[l + 4a/3V] 
a In o" J 



e~ 2a/3y °Erfc 



2a : 



■(^o-/3a 2 ) 



a 3 / 2 a 



(3.24) 



U 3 (a) [A^f3a + 8 (a) 3 / 2 (f3a) 3 - 4a 3/2 /3 2 aF ] 



+ 



dU 3 (a) 
dlna 



[2^a~f3a + 2a 3/2 (/3a) 3 - 2a 3/2 f3 2 aY ] 



As in the case of the leading order term we notice the presence of additional mixing 
terms in powers of (3. This clearly indicates that contrary to standard derivations of 
the NG mass function based on the Press-Schecther formalism, the effect of the non- 
spherical collapse of halos cannot be simply reabsorbed in the Gaussian part of the 
mass function, but directly alters the non-Gaussian dependence 6 . 



4 Local and Equilateral Non-Gaussian Mass Function 
4.1 Reduced bispectra and fitting functions 

In deriving the non-Gaussian terms we have made no assumptions on the specific type 
of non-Gaussianity. In the case of non-Gaussianity with non-vanishing skewness the 



6 Notice that the non-Gaussian terms for (3 = do not exactly coincide with the corresponding 
formulae in [50] neglecting the mixing terms in an and Jnl- The mass function for a diffusive 
stochastic barrier cannot be inferred from the spherical collapse scale by a simple rescaling of the 
variables including the skewness of the initial density field. 
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three-point correlation function of the smoothed linear density field at z = formally 

reads as: 



x M(A; 1 )M(A; 2 )M(A;3)(C(k 1 )C(k 2 )C(k3)) c , (4.1) 

where 

Mk) = —^-T(h)k\ (4.2) 

with Hq being the Hubble constant, Q m the matter density relative to the critical one, 
T(k) the transfer function at z = and ((k) is the curvature perturbation with 

(C(k 1 )C(k 2 )C(k 3 )) c = (2n) 3 5 D (k 1 + k 2 + k 3 )B ( (k u fc 2 , fc 3 ), (4.3) 

where B^(ki, k 2 , k 3 ) is the reduced bispectrum and encodes the characteristic momen- 
tum dependence of the non-Gaussian signal. The dependence on the Dirac-5 function, 
5o(k), ensures the momentum conservation, implying that admitted momentum space 
configurations form closed triangles. 

A variety of models, including standard inflation, predicts 

Bdh: k 2 , h) = [P c (fc 1 )P c (fc 2 ) + P ( (h)Pdh) + P c (k 2 )P c (h)} , (4.4) 

where f l ^ c L is a constant amplitude and P^{k) = A s k ns ~ A is the primordial power 
spectrum. This type of non-Gaussianity is also known as "local" since the curvature 
perturbation at any point in space can be written in terms of a linear Gaussian field, 
Cg(x): 

C(x) = Cg(x) + [C G (x) - (C G (x))] . (4.5) 

In such a case the non-Gaussian signal peaks around momentum space configurations 
for which one of the momenta is much smaller than the others (e.g. k\ <C k 2 ~ fc 3 ). 
WMAP-7yr measurements of the CMB bispectrum constrain the amplitude of local 
non-Gaussianity in the range — 10 < f^Q < 74 at 2a [58]. 

Alternatively, models in which primordial density fluctuations are generated by a 
field subject to higher-derivative terms predict a reduced bispectrum of the form 



B c (k u k 2 ,k 3 ) = ^f N tA 2 



+ 



5 jnl s[ 2 kt na kt na 3(fc 1 fc 2 A;3) 2 ( 4 -™»)/ 3 
5 perm. 



(4.6) 



This is also known as "equilateral" non-Gaussianity, since the signal maximizes in 
momentum space configurations for which ki ~ k 2 ~ k 3 . WMAP-7yr measurements 
constrain the amplitude of equilateral non-Gaussianity in the range —214 < f^ 1 < 266 
at 2a [58]. 
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Figure 2. £3 (a) and U3 (b) for local non-Gaussianity. Black solid curves are the numeri- 
cally computed functions obtained from numerical integration of eq. (4.1). Red solid curves 
correspond to fitting formulae eqs. (4. 7), (4. 9) with a few percent errors (bottom panels). 




S = a z S = <jZ 



Figure 3. As in figure 2 for equilateral non-Gaussianity. Fitting formulae are given by 
eqs. (4.11)-(4.13) 

We integrate eq. (4.1) numerically for local and equilateral type of non-Gaussianities 
and infer £3, dS^/d In a, U3 and dlA^/dhia. The transfer function has been computed 
for a vanilla ACDM model best fit to WMAP-7 years data using the CAMB code [64]. 
The functions are plot in figure 2 and 3 respectively, over the mass range probed by 
[38], which corresponds to —0.2 < ln(l/cx) < 0.8. These functions can be approximated 
to better than a few percent accuracy by simple fitting formulae. In the local case we 
find 
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1 5fi 

S l r(S) = /iVL^lO" 4 (4.7) 
dSt _ 6.314985 

- £1.4736 iU 

W 3° C ( 5 ) = /^^S 10 " 4 ( 4 - 9 ) 
d^ oc „ 1.0122 



dS J S 1A 
while for equilateral NG we have 



/tvl-^10" 4 (4.10) 



ST\S) = -f NL 5 -§^lO- 5 (4.11) 



dS c 3 qui „ 2.5 



Wf ui (S) = -/^^10" 5 (4.13) 



dW 3 equi „ 3.108 



4.2 Validity of the perturbative expansions: limits on f NL 

We have inferred the halo mass function in the presence of primordial non-Gaussianity 
as an expansion about the Gaussian Markovian solution by Taylor expanding the path- 
integral eq. (2.13) to first order in the three-point correlation function. In this section 
we derive limits on Jnl that bound the validity of such a perturbative approach by 
imposing the ratio |/vg(<t)//o(«t)| << 1 with f NG (a) = fk G ( a ) + /jvgO 7 )- 

In figure 4 we plot /tvg( ct )//o( ') f° r local (left panel) and equilateral (right panel) 
non-Gaussianities in units of /jvl for the barrier model parameters shown in figure 1. 
First, we may notice that in the local case the ratio increases toward large value of 
log 10 a -1 and maximizes in the high-mass end. Secondly, the overall amplitude of this 
ratio varies for different barrier model parameter values. This is an indication that 
the non-Gaussian imprint on the halo mass function depends on the collapse model 
which we will discuss in more detail in the next section. A similar trend occurs in the 
equilateral case, the main difference result in the low-mass behavior and the value of 
the ratio in the high-mass end. 

From these curves we infer the ratio /jvg(c)//o(c) ~ 0.5 for \f^1\ ~ 1000 i n the 
case of the Diffusive Barrier while for rigid barriers, (i.e. Db = 0), we find ~ 
170. For the equilateral type the limits are l/j^l ~ 1000 for the rigid barrier and 
l/ivxl ~ 2500 for the diffusing barriers respectively. 
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Figure 4. /atg(°')/./'o( ') as a function of log 10 ((j _1 ) for local (top panel) and equilateral 
(bottom panel) non-Gaussianities in units of /jvx- Different curves correspond to the collapse 
parameters values of models shown in figure 1 



The mass function derived here can be used to compared with N-body simulation 
and for observational analysis purposes to test primordial NG for a large amount of 
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primordial non-Gaussianities and mass range. 

5 Halo Collapse Model and NG Signal 

We now discuss the imprints of primordial non-Gaussianity on the halo mass func- 
tion and the dependence on the halo collapse model parameters. In order to evalu- 
ate the NG features it has become standard to consider the ratio 1Z = /(cr)//c(cr), 
which we plot in figure 5 for local NG (top panel) and equilateral NG (bottom panel) 
with fisri = —80,-50,50,80 in the case of non rigid barriers (Db = 0) and /atl = 
—800, —80, 80, 800 in case of diffusive ones. We can see that independently of the 
value of (3 and Db the deviation from the Gaussian case maximizes in the high-mass 
end consistently with expectations that NG affects primarily the large masses both 
in the local and equilateral case. However, we also notice that the non-Gaussian sig- 
nature depends on the barrier model parameter values. In particular, by comparing 
curves with the same value of Jnl and different values of j3 and Db we can infer two 
distinctive trends. On the one hand the average drift of the barrier enhances the NG 
imprint, especially at small masses. On the other hand, the barrier diffusion strongly 
supress the overall amplitude of the non-Gaussian signal relative to the Gaussian one, 
particularly in the high-mass end. Similar trends occur in the equilateral case. 

As already inferred from the form of the non-Gaussian terms in section 3, this 
shows that the effects of the statistics of the primordial matter density field and the 
non-spherical collapse of halos are coupled. 

Furthermore the nature of the imprints is consistent with the excursion set picture. 
In fact, primordial non-Gaussianity manifests as an additional non-Markovian correla- 
tion to the random walk trajectories. The barrier diffusion increases the likelihood of 
first-up crossing, while the average drift reduces it and as discussed in section 2.4, the 
former has the effects of diminishing the amplitude of the non-Markovian term relative 
to the Markovian one, while the latter enhances it. 

Furthermore, the fact that the non-Gaussian signal depends on the halo collapse is 
particularly important from an observational point of view. In fact, if we don't use the 
right halo mass function we might put wrong prior on the amount of non-Gaussianity. 
It is plausible that such degeneracies may be broken including measurements of the 
redshift evolution the halo mass function. Although, this will require predicting the 
redshift and cosmology of the barrier model parameters for non-Gaussian initial con- 
ditions, a study which goes beyond the scope of this paper. 

6 NG Halo Mass Function and N-body Simulations 

In recent years, several groups have performed numerical N-body simulations with 
non-Gaussian initial conditions to probe the imprints of primordial non-Gaussianity 
on the distribution of dark matter halos [36-39]. These studies have provided fitting 
functions for the non-Gaussian halo mass function and tested the validity of analytical 

7 For instance we performed the analysis for 2.10 13 < M < 2.10 15 /i _1 M Q at z = 0. 
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Figure 5. TZ = /(o")//g(c) as function of log 10 (<r~ 1 ) for local (top panel) and equilat- 
eral (bottom panel) non-Gaussianities using different values of the collapse parameters and 
different non-Gaussian amplitudes. 

formula based on the Press-Schechter formalism. Here, we confront the mass function 
derived in the previous sections against the results of Pillepich, Porciani and Hahn [38]. 
These authors have performed a series of high-resolution runs with 1024 3 particles 
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in large volume boxes resolving halos with mass 2 x 10 13 < M[h M Q ] < 2 x 10 15 
at z=0 for local NG initial conditions in the case of WMAP-3yr and WMAP5-yr 
calibrated cosmologies. Halo mass functions have been measured using the Friend-of- 
Friend algorithm with linking length b = 0.2 at different redshifts and fitted with the 
standard parametrization 



in doing so the authors have inferred polynomial fits which parameterize the depen- 
dence on Jnl and z of the coefficients A, B, C and D. Here, we limit our theoretical 
model comparison with the results at z — 0. In what follow we use the parameters of 
Table 5 of [38] for which the accuraty is of order 5% on the range —0.2 < h^cx" 1 ) < 0.8 

Before discussing the results of this comparison let us make a few remarks. The 
excursion set mass function for the diffusive barrier model depends on two input param- 
eters which specifies the statistical properties of the stochastic barrier. Other quantities 
such as the spherical collapse threshold S c and the amplitude of the filter-correction 
k are predicted for a given cosmological model. For instance in ACDM WMAP5-yr 
cosmology 5 C = 1.673 at z = and k = 0.475. 

In principle, the values of /3 and Db can be predicted for a given cosmology. As 
already stressed in section 2.2 this requires estimating the properties of the stochastic 
barrier by numerically solving the ellipsoidal collapse model equations. As shown by 
Doroshkevich [60], the randomness of the initial density field causes the parameters 
describing a homogeneous ellipsoid to be random variable themselves with a character- 
istic probability distribution. By generalizing this to non-Gaussian initial conditions 
and solving the ellipsoidal collapse model equations one can infer the probability dis- 
tribution of the collapse threshold as function of the variance of the linear density 
field. Then, the stochastic barrier model parameters can be derived by matching their 
value to the moments of the probability distribution. Such an inference will provide 
a complete prediction of the cosmology, redshift and non-Gaussian dependence of (5 
and Db which at moment is missing. However, testing such dependences against nu- 
merically inferred N-body halo mass functions will first require addressing the effect of 
the second order cloud-in-cloud problem, which as we have mentioned in Section 2.2 
manifests as an additional scatter on the halo collapse threshold. 

In the absence of theoretical prediction for the values of (3 and Db, the diffusive 
drifting barrier model remains an effective description of the non-spherical collapse of 
halos. Consequently, the model parameters can be determined only through compar- 
ison with N-body simulation results. For instance, in [51, 52] we have confronted the 
Gaussian mass function N-body simulation data by Tinker et al. [63] obtained using 
the Spherical Overdensity (SO) algorithm. This has allowed us to infer the best fit 
values Pq° = 0.057 and Dg = 0.294 and quite remarkably found that for such pa- 
rameter values the mass function reproduces the data to numerical accuracy (~ 5%). 
However, we cannot assume these values to confront with the results from the Pillepich 
et al. simulations even for f^i = 0. This is because there is roughly a ~ 10% scatter 
between SO and FoF mass functions due to systematic differences of the halo detection 
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Figure 6. Relative difference of the NG mass function with respect to Pillepich et al. fitting 
formula for the Gaussian calibrated barrier model parameters. The discrepancy remains 
within the numerical accuracy of the N-body results (solid parallel lines) for — 10 < /jvl ^ 30. 

algorithms. Therefore we have calibrated the value of (3 and Db against the Gaussian 
mass function from [38]. 

We find /3 G = 0.0529 and Dg = 0.328. Assuming such values we plot in figure 6 
the relative difference between eq. (2.19) and eq. (6.1) for different values of f NL . We 
can see that for Jnl = the differences are well within the numerical uncertainties 
~ 5%, which is consistent with the results of the comparison to Tinker et al. [63] 
presented in [51, 52]. As it can be seen in figure 6, for ^ the relative difference 
remains within the numerical errors only in the range —50 < /jvl < 150, while for 
larger values the mass function evaluated using the Gaussian calibrated parameters 
largely deviates from the numerical fit. 

This simply implies that for J'nl outside this range the non-spherical collapse of 
halos strongly differ from the Gaussian case and consequently we can expect (3 and D b 
to vary with f^L- 

We can see this more clearly by inferring the best fit values of f3 and D b for differ- 
ent values of /nl in the range —80 < /jvl < 750 for which we compare equation (2.19) 
against eq. (6.1). We may notice that (3 &t increases for increasing values of /jvl > 0, 
as well as Dg. 

In particular, for large positive values of f^L the average drift is about 4 times 
larger than in the Gaussian case. This indicates that for large primordial non-Gaussianity 
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Figure 7. Best fit values of j3 (top panel) and Db (bottom panel) relative to the Gaussian 
case as function of /wi- 



the collapse of dark matter halos is highly non-spherical. Furthermore, since the drift 
coefficient increases for increasing values of Jnl, halos with higher mass will be privi- 
leged relative to the smaller ones. 

Similarly, the diffusion increases with respect to the Gaussian case for increasing 
values of /jvx- It is worth reminding that the barrier has a Gaussian distribution of 
mean value B = 5 C + j3S and variance Db- Therefore, for positive value of Jnl, the 
non-Gaussian nature of the initial density field tends to erase the imprint of the mean 
ellipsoidal collapse value which carries the signature of the gravitational dynamics. 

It would be interesting to perform a study similar to that of Robertson et al. 
[65] to infer the distribution of the linear collapse threshold of halos detected in non- 
Gaussian simulations, which can provide a direct confirmation of the trend of /3 and 
D B as function of /jvx obtained here. 



In figure 8 we show the relative difference with respect to Pillepich et al. fitting 
formula of the NG mass function computed with the best fit barrier parameters (red 
long-dash dot line) for the values of f l fi c L used in figure 7. We can see that the differences 
remain within the numerical accuracy of the N-body results. In addition, we show the 
relative difference for the Loverde et al. [43] formula (LV, green dot line), Matarrese- 
Verde-Jimmenez [42] (MVJ, magenta short-dash line) and Lam-Sheth [44] (LS, cyan 
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Figure 8. Relative difference with respect to the Pillepich et al. fitting formulae for different 
values of f^§ c L of non-Gaussian halo mass functions for the DDB case with best fitting barrier 
parameter values (red long-dash dot line), Loverde et al. [43] (green dot line), Matarrese- 
Jimenez- Verde [42] (magenta short-dash line) and Lam-Sheth [44] (cyan short-dash dot line) 
respectively. 

short-dash dot line). Since the Gaussian limit of these different formulae is given 
by the Press-Schechter mass function, for a comparison with the N-body simulation 
results we have renormalized them to the Gaussian fitting formula by Pillepich et al. 
[38]. We find these formulae to be close to the numerical accuracy of the N-body 
results for —80 < f l £ c L < 150. As we infered in from 5 this correspond to the range 
where the non-Gaussian dependence of the non-spherical collapse is negligible. In 
the light of this analysis we find that our halo mass function explains the so called 
"fudge" factor which is used to fit the non-Gaussian simulation. Moreover outside 
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the range -50 < /#£ < 150 the non-spherical nature of the halo collapse must be 
taken into account when testing primordial non-Gaussianity with halo mass function 
measurements. 

7 Conclusion 

In the upcoming years several observational program will provide precise measurements 
of the number counts and spatial distribution of galaxy clusters. These measurements 
will provide new insights on the mass distribution of dark matter halos and conse- 
quently on the underlying cosmology and the statistics of the primordial density field. 

Here, we have performed an ab initio calculation of the NG halo mass function 
using the path-integral formulation of the excursion set formalism. We have considered 
a stochastic model of the halo collapse threshold which captures the main features of 
the non-spherical collapse of halos. We have specifically focused on primordial NG 
due to a non-vanishing skewness of the initial density field and inferred the mass 
function through a perturbative expansion of the path-integral. The computation can 
be extended to higher-order correlations and we have showed that the effect of a non- 
vanishing kurtosis can be easily obtained at leading order. 

Contrary to standard approaches based on NG extension of the Press-Schechter 
calculation, we find that the effects of the non-spherical nature of the collapse of halos 
on the mass function directly couples to the signature of primordial non-Gaussianitiy. 
As such, these effects cannot be simply reabsorbed in a Gaussian normalization term. 
We have studied the case for local and equilateral non-Gaussianities, though the for- 
mulae provided here can be applied to any NG shape function. 

We have compared the inferred mass function with NG local N-body simulations 
results and found a remarkable agreement. Using the barrier model parameters cali- 
brated on the Gaussian N-body mass function we find that for large non-Gaussianities 
(/jvi < —50 and /jvl > 150) the non-spherical collapse largely deviates from the Gaus- 
sian case and its effect cannot be simply reabsorbed in a Gaussian normalization as in 
the case of PS derived mass functions. These effects can, therefore, be important for 
cluster count studies since the non-Gaussian dependence of the non-spherical collapse 
of halos changes the imprint of primordial non-Gaussianities on the mass function. For 
instance in [66] the author has shown that the abundance of clusters which do not be- 
long to the super clusters, detected in Gaussian N-body simulations, is well described 
by the Gaussian mass function derived in [51, 52] for Db = and a non- vanishing j3. 
In light of this result we can expect that the number counts of this type of clusters 
may be more sensitive to the signature of primordial NG since the imprint has a larger 
amplitude compared to field clusters with D B ^ and (3 = 0. 

These results suggest several directions of future investigation. On the theoretical 
side it appears evident that a theoretical prediction of (3 and Db based on a detailed 
study of the ellipsoidal collapse model is needed. This will allow us to predict the cos- 
mology and redshift dependence as well as their relation with respect to statistics of the 
primordial density field. Thus, a barrier model parameter inference using mass func- 
tion measurements may provide a way of distinguishing between different cosmological 
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scenarios assuming a proper treatment of the second order cloud-in-cloud problem. On 
the numerical side it will be interesting to confront these model predictions with the 
linearly extrapolated collapse threshold of halos detected in N-body simulations for 
different cosmological scenarios. 
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